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Abstract 

The idea of generating integrals analogous to generating functions is first 
introduced in this paper. A new proof of the well-known Finite Harmonic Series 
Theorem in Analysis and Analytical Number Theory is then obtained by the 
method of Generating Integrals and Fractional Calculus. A generalization of 
the Riemann zeta function up to non-integer order is derived. 



1 The Finite Harmonic Series Theorem 

Definition 1 The finite harmonic series is 

n 1 1 1 

fc=l 

Theorem 1 (The Finite Harmonic Series Theorem) 

It is well known p. 16, (1.7.9)] that 

h{n) = V(l+™)+7 (1) 
= logn + 7 + 0(l/n) (2) 

d T'(z) 
where the digamma function i^iz) = — log(r(z)) = — — -, the Euler constant 

dz T(z) 

7 = — r'(l) = 0.577215 and the prime denotes differentiation. 

The first part (|l]) of the Theorem can be proved by differentiating the recur- 
rence relation of r(n+l) = nT(n) and summing over n |Q p. 256]. The second 
part (|2|) can be proved by using the Euler-Maclaurin summation formula [Q, p. 
629]. 

In this paper, we shall present an alternative elementary proof |J of the 
first part ([!]) of the Theorem by using Fractional Calculus and the idea of 
generating integrals. 
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2 Introduction Sz Motivation of Fractional Calculus 



"Thus it follows that d}l 2 x will be equal to . . .from which one day 
useful consequences will be drawn. " - Leibniz in a letter |5J to 
L'Hospital. 



"When n is an integer, the ratio d n p, p a function of x, to dx n can 
always be expressed algebraically. Now it is asked: what kind of ratio 
can be made if n is a fraction?" — Euler ||. 



"The idea of an integral or derivative, of arbitrary non-integral or- 
der, was introduced into analysis by Liouville and Riemann. Such 
integrals and derivatives may be, and have been by different writers, 
defined in a variety of manners, and different systems of definitions 
may be the most useful in different fields of analysis. " — Hardy and 
Littlewood 0. 



In this paper, we shall consider D° , a € R, an operator of non-integer 
order, a notion first pondered upon by Leibniz |, pp. 301-302], Euler g p. 55], 
Lagrange ||, Laplace |], p. 85 and p. 186], Fourier |nj and Abel Jll| during the 
late 17th century to early 19th century. We shall briefly review several known 
and equally valid definitions for D a , and then focus on one of the definitions, the 
Riemann-Liouville (R-L) Fractional Calculus. The foundation of R-L Fractional 
Calculus was laid by Riemann [jT2| and Liouville |l3| in the late 19th century, 
and then subsequently developed by Cayley [fH|], Laurent |j~5[ |, Heaviside [16|, 



Hardy and Littlewood |7L 17, 18, M, and many others. It was largely regarded 



as a mathematical curiosity until only recently when Mandelbrot, the discoverer 
of Fractals, found an application of the R-L Fractional Calculus in the Brownian 
motion in a fractal medium, and speculated a possible connection between the 
analysis of Fractional Calculus and the geometry of Fractals [2C]. 

We shall also introduce the idea of generating integrals by analogy to gener- 
ating functions. As we shall see, just as a certain generating function is useful 
for generating a certain desired sequence of numbers, a certain generating inte- 
gral is similarly useful. While a generating function f{z) generates a sequence 
of numbers {p n } in the coefficients of the terms of different orders in its power 
series expansion ^ n p n z n , a generating integral generates a sequence of numbers 
q n in the coefficient of a term in the result of an n-fold integration. 

However, a generating integral has one unique advantage over a generating 
function. The R-L Fractional Calculus can be used to analytically extend a 
generating integral of iteration order n € Z + to order p € R. The result is that 
the sequence of numbers {q n } is in turn analytically extended to a function 
q(p), pGR. 
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We shall then show how the Riemann-Liouville Fractional Calculus and 
the idea of generating integrals can be used to prove the well-known Finite 
Harmonic Series Theorem. 



3 Differential-Integral Operator D r 



Definition 2 Let the operator D^ a , n G Z, acting on a function f at the point 
x be defined as 



DXf(x) = 



in 



/(*) 



(n > 0) 
(n = 0) 



(3) 



f X f(x)(dx)- n (n<0) 

V Ja 



where the n-fold integration is defined inductively as 



rx rx n «„-i rxz rX2 

f(x)(dx)~ n = / / •••/ / f{xi) dxi dx 2 ■ ■ ■ dxn-! dx n . 

a J a J a J a J a J a 



n-times 

As an example, consider f(x) = x r ' 

r o 



x m~n = T(l + m) 

(m—n)\ r(l+m — n) 



, un (\m\-l+n)\ n 
1 j (H-l)! 

r(l+m+e) 

= hm — ; 

e^o T(l+m+e — n 



X 

a 

~ rt 
■ X 



(m > 0, m < n) 
(m > 0, m > n) 

(m < 0, m < n) 



la [la ^ ^) ^ 



(m—n) 



(m < 0, m > n) 



= (-1) 



m+l 



N-l)! 7a 



logx((ix) < - 



m—n) 



where x r ' 



(4) 



= (x m " n -a m ~ n ). 



If we tabulate Z7| a x m for n, m G Z and omit the constant terms containing 
a, we can observe a pattern emerges as in Table 1. 
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L»"x 


m 
















-3 


-2 


-1 





1 


2 


3 


2 


2!/5! x 5 


2!/4! x 4 


2!/3! x 3 


x 2 


2! 


se 


2! 









1 


1/4! x 4 


1/3! x 3 


1/2! x 2 


X 




1 















1/3! 


1/2! x 2 


X 


1 















-1 


/logx (dx) 2 


/ log x (dx) 


log X 


x _1 


- X 


-2 


2! x" 3 


-3! 


X~ 


-4 


-2 


- / log x (dx) 


- log X 


- x _1 


x- 2 


-2! x" 


-3 


3! x" 4 


-4! 


af 


-5 


-3 


1/2! logx 


1/2! x" 1 


-1/2! x" 2 


X -3 


-3!/2! ar- 


-4 


4!/2! x" 5 


-5!/2! 


X~ 


-6 



Table 1: Tabulated results of D™x m 

4 Riemann-Liouville (R-L) Fractional Calculus 

The R-L Fractional Calculus pl| begins with 



/x rx„ rx n -i [-X3 rX2 
/ / "'/ / f(x\) dxidx 2 - ■ ■ dx n -idx n (5) 
Ja Ja Ja Ja 



n-times 



for n E Z+ as the fundamental defining expression, and it can be shown [21, p. 
38] to be equal to the Cauchy formula for repeated integration, 



1 



/(*) 



T(n) Ja (x-t) 



1-71 



dt . 



(6) 



Definition 3 (R-L Fractional Calculus) 

The R-L fractional integral is analytically extended from Ij^j as 



d a 



dx 



f(x)(dx) CT by extending (3) 



1 r* f(t) 



T(-a) J a (x-ty+ 



dt (a < 0, a, a G R) by 



(7) 



and the R-L fractional derivative is in turn derived from the R-L fractional 
integral (jj) by ordinary differentiation: 



D: {a f(x) = D™ a (D^- a) f(x)) (a > 0, m G T 
where m is chosen such that m > 1 + a, a > 0. 



(8) 



Lemma 1 The equation ^) is independent of the choice of m for m > 1 + 
a, m G Z+, a G R, a > 0. 
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Proof 

For m > 1-Hr, m € Z + , a > 0, we have — (to— a) < — 1 < and (m — <r — 1) > 
0. The first condition, — (to — a) < 0, allows us to use the equation (^) to write 

i r f(t) 



D- {m - a) f(x) 

x a J \ / 



r(m-o-) A (x-t) 1+,T - m 



From ® 



1 



/(*) 



0?X m Vr(TO-Cr) J a (x-t) 1+c7 - m 

1 P ( d m 



dt 



r(m-cr) 7a 



/(*) 



■(x-t) 



m— a— 1 



dt 



0) 



The second condition, (to— a— 1) > 0, and the condition to > allow us to use 
the second case of @. Thus, (||) becomes 



r(m-cr) Ja 
1 



/(*) 



T(to — a) 



T(-a) 
x fit) 
r(-a) J a (x-ty+° 



(x-t)~ {l+ ^ dt 
dt 



□ 



When a = 0, (]?]) for f[x) = x r is well-defined only for the half plane r > — 1. 
Consequently, in the R-L Fractional Calculus, D^ a x r is well-defined only for 
the half plane r > — 1. 



T(l+r) 
r(l+r-<r) 

T(l + r) 

L r(i+r-ff) 



(<t > 0, r > -1) 

(<t = 0, Vr) 
(<7 < 0, r > -1) 



(10) 



5 Fractional Calculus by Cauchy Integral 



The Cauchy Integral for an analytic function f(z) in the complex plane [22, p. 
120] is 



f (n) (z ) 



T(l + n) 
2m 



c {z- z ) 1+ ™ 



dz 



(11) 



Analytic extension of the Cauchy Integral from n E Z + to s £ R + gives an 
analytic extension of D n f(zo). However, the analytic extension is not trivial. 
The term (z — zq) 1+(T will become multi- valued and thus the result may depend 
on the choice of branch cut and integration path. 
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6 Fractional Calculus by Fourier Transform 



In the theory of Fourier Transforms, 

/>) = / + °° f(x)e lkx dx, 



— oo 

oo 



fix) = ±- I "' f(x)e-* kx dk, 



2tt 



and 



oo 
-oo 



L%f{x) = / + °° f(k)I%(e- ikx )dk (aeR) 

{-ikff{k) e~ ikx dk . 
This approach is often known as the pseudo-differential operator approach. It 



was shown by Zavada [23] to be equivalent to the Riemann-Liouville Fractional 



Calculus and the Fractional Calculus by Cauchy Integral. 



7 Functional Analytic Approach 

In the functional analytic approach, an example of a functional integral of an 
operator A is 

= _ fV" 1 (Al - A)' 1 A dX (0<o<l, AgR). (12) 

vr Jo 

(Al — A) is called the kernel of the functional integral. The evaluation of 
the integral with respect to real variable A requires various conditions on the 
spectrum of the operator A. 

The analytic extension of D in the functional approach is then obtained 



from replacing A by D in (12). 



For details of this well-developed functional analysis approach, see [24 



8 Differentiating and Integrating in non- integer s- 
dimensions 

The differential of an integer n-dimensional function in n-dimensions can be 
expressed as 

d d d 

a — f{xi,x 2 ,---,x n ) . 

OX\ OX2 ox n 

The corresponding integral can be expressed as 

f(xi,x 2 ,---,x n )d n x= •••/ / f(x!,x 2 , ■ ■ ■ ,x n ) dx\dx 2 ■ ■ ■ dx n . 



n-times 
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If / is spherically symmetric, / = f(r), then 

d d 8 d 11 ' 1 d _ r(n/2) d n ~ l 



dx 1 dx 2 dx n J dr™" 1 dQn-i 2vr n / 2 dr n - lJ v ' ^ 
can then be analytically extended to the differential of a non-integer s-dimensional 
function in s-dimensions, 



&S ~-- f(r) = r(5/2) ^ f(r) 



dr 3 - 1 dVls-i ' 2tt s / 2 dr s 
where s € R or s € C, and f2 is the n-dimensional solid angle. 
Similarly, the corresponding integral 



(14) 



/ d n x 



oo r2ir 

r n - l f{r)dr J dO x / sin0 2 d02 • 



sin n 2 6» n -id6> n _i 



o Jo 
r n - l f(r)dr / d^-l = — — / r n ~ l f{r) dr 



can be analytically extended to 



2tt s / 2 f°° 
f d s x = — — — / r s ~ x f(r)dr . 

T(s/2) y v ; 



(15) 



(16) 



This method was developed by 't Hooft and Veltman [25| in 1960's. The 
method was central to an important technique called Dimensional Regulariza- 
tion in Quantum Field Theory where the method is used to isolate singularities 
in divergent integrals. 



9 Generating Integral of the Finite Harmonic Series 

Theorem 2 (Generating Integral of h(n)) 



h(n) = log x 



r(l+n) f x 



[ \ogx{dx) n (neZ) 
jo 



(17) 



where 



m (dx) n 



x rx n fX n -i rx :j px 2 



a Ja, Ja Ja Ja 



f(x\) dx\ dx2 ■ ■ ■ dx n -i dx r 



n -times 



Proof 

We observe that —h(n)/n\ appears in the coefficient of the x n term when 
we repeatedly integrate logx: 



/ log x (dx) = x( log x — 1 ) , 
Jo 

[ logx(dx) 2 = %-(logx- -) , 
Jo 11 



\ogx(dx) n = — -(log x — h(n)) 



(18) 
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We can prove this observation by induction: 



log x (dx 



,n+l 



o wo 



\ogx(dx) n (dx) 



-(log x — h(n) ){dx) 
o n\ 

-Alogx-hin)) - I ——-(dx) 



(n + l)\ 



„n+l 



(n + l)\ 



log x — h(n) 



o (n+1)! 
1 



n + 1 



(logx — h(n+l) ) . 



(n+1)! 

Rearrangement of (|i"8|) yields the Theorem. 



□ 



Theorem 3 (Generating Integral of h(p)) 



h(p) = log x 



r(i+p) 



/ log x (dx) p (p G R) 



(19) 



Proof 

By analogy to generating functions, we take 

log x (dx) n 

as the generating integral of the finite harmonic series h(n), and so the natural 
analytic extension of the generating integral takes the form of 



log x (dx) p = — 

o y 1 r(i+p) 



(logx - h(p)) 



Noting that log x may be expressed as 

rx rx ^ 

logx = / x^dx = lim / x~ l+e dx = lim - (x e — 1) 
J l e->0 J 1 ^0£ l ' 



(20) 



(21) 



we can now evaluate the fractional integral in ( |20| ) by the R-L Fractional 
Calculus (Definition ||) . 

/ logx (dxY = lim - (x e - 1) (dx) p = lim - / (x £ - 1) (dx) p 
Jo Jo e e-*o e Jo 



lim — 

e^O e 

lim — 

e^O e 

x p 
lim — 



x e (dx) p - / l(dx)' 5 
o Jo 



r(l+e) x e 



1 



r(i+e+ P ) r(i+ P ) 



(22) 
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where the interchange of the integral and the limit is justified by Arzela's the- 
orem on bounded convergence [26, pp. 405-406] as (x e — l)/e is integrable in 
x € [0, x] . 

Combining p3) with ( 20 ) gives the analytic extension of the finite harmonic 



series: 



Hp) 



log x — lim — 

£-0t 



r(i+e)r(i+ P ) 



x e - 1 



log x — lim — (x e 



+ lim — 



1) 



r(i+e)r(i+ P ) 



T(l+e+p) 

r(i+6)r(i+ P ) - 
r(i+e+p) . 



log X 



lim — 

e— o e 



r'(i+p) 



lim — (x e 



+ lim — 



r(l+e)r(l+p) 
r(l+e+p) 



r(i+e)r(i+ P ) 
r(i+e+p) 

-r'(i) 



(23) 
□ 



" r(i+ P ) 

= V(l + P) +7 
where the limit has been taken with L'Hospital rule. 

3r 

2.5 
2 
1.5 

%) i 
0.5 




Figure 1: The curve h(p) = + p)+7 passes through the points (n,h(n)) 

n 1 

where n € Z + , h(n) = ^ 




k=l 



k ' 



Hence, by the application of the R-L Fractional Calculus to analytically 
extend the generating integral, we have found an alternative elementary proof 
of the first part (||) of the Finite Harmonic Series Theorem. 



Theorem 4 



log x (dx 



t(i+p) 



(logx- V(i+/?) -7) (peR) 



(24) 
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Proof 

Replacing the h{p) in (^) by ( p3| ) gives the Theorem. 



□ 



10 The Riemann Zeta Function up to Order n 

The analytic extension (^3|) can be generalized. 

Definition 4 The Riemann zeta function 

1 

k 



CW = E tJ ( Re ( s ) > 1)» (») 
fc=i 



i/ie poly gamma functions p. 260, (6.4-1)] 

(26) 

anc? i/ie Riemann zeta function up to order n 

n -1 00 -. 

C(- I ") = E TJ = CW - E p ( Re ( s ) > !) ( 27 ) 

fc=l ^ fc=n+l * 

may 6e combined to write 

C(™ I *) = Epr = (b^Tyr (y-^a+n) - ^ {m - 1} (i)) 

f-lf d m x=n 

= ^-—log(T(l+x)) (28) 
(m — lj! ax m x=o 

The analytic extension is then obtained by replacing the derivative in (ffif) 
with a fractional derivative: 



C(*|*) = ^ 2£ log(T(l+20) * * (.s.ceC) (20) 

1 (s) z=0 



which can be evaluated when log(r(l+x)) is expressed in the form of an asymp- 
totic series [||, p. 257, (6.1.41)]. However, w(s) depends on the choice of exten- 
sion to the R-L Fractional Calculus (Definition ||) into the other half plane, 
r < -1, cr € R. 



11 Analytic Extension of R-L Fractional Calculus 

Consider the case of D™ where a = 0. We shall introduce the (cr, r) diagram in 
which the numerical factor of D£x r is mapped to the point at coordinate (cr, r) 
of the diagram. The (cr, r) diagram of D£x r can be characterised into 4 regions 
as in Figure ^: 
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Definition 5 (Regions of (cr, r)) 



The zero region 


Zero 


= {(*, 


r) 


: r < a, 


r 


>o} 


the upper region 


Upp 


= {(*, 


r) 


: r > a, 


r 


>o} 


the lower region 


■Low 


= {(*> 


r) 


: r < a, 


r 


<0} 


the log region 


Log 


= {(*■ 


r) 


: r > a, 


r 


<0} 



A point lying on the right of the r-axis (cr > 0) is a differentiation; a point 
on the left (a < 0) is an integration. 

The (cr, r) diagram at integer grid points everywhere except in the log region 
gives numerical factors identical to those in Table 1. 

An extension of R-L Fractional Calculus to the other half plane r < — 1 is 
given by 

logx in f2 



e^o T(e) x 



lim 



r(l+r+e) 



(30) 



elsewhere 



e->o T(l+r+e-a) 
{(cr, r) : cr € R, r £ Z~}, the set of horizontal lines in lower and log 



where fi 
regions. 

r(l + r)/r(l + r — a) is finite everywhere in the zero and upper regions. 

e — > r(l + r + e)/r(l + r + e — cr) is well-defined everywhere in the lower and 
log regions except in f2\(Z~ x Z"). Following Theorem |4], the R-L fractional 
integral of logx can be evaluated exactly and expressed in only elementary 
functions, 



x 



r(l-a) 

T(l+r+e 



(log a; --0(1-0-) -7) (cr € R\Z) 

(cr € Z) 



(31) 



lim 

e^o r(l+r+e-cr 
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(j3C|) is thus well-defined. 

To analytically extend from D°x r on the real plane (a, r) to D£x r on the 
product of complex plane and real line, (s, r) € C x R, we simply replace a € R 
in @) by s G C. 

12 Open Problems 

1. Generalize ( j30[ ) for D^_ c z w , s,w,z,c € C. For the complex function z™, 
one has to specify, in addition, the integration contour for Re(s) < 0. 

2. Find the exact expression for log(r(l+x)) and w(s) in terms of elementary 
functions in analytic extension of the R-L Fractional Calculus given by 

(ED- 

13 Tables of Generating Integrals 

"Nature laughs at the difficulties of integration. " — Laplace 

An interesting consequence is that objects of the form J x r (Jog x) a (dx) p , 
r,a,p G R, exist and can be generating integrals for certain functions. Perhaps 
it may be worthwhile to introduce, in the future editions of Tables of Integrals, 
a new section which gives the coefficient of the x k term, w(p,r,a,k), corre- 
sponding to these generating integrals to facilitate the evaluation of integrals 
of similar forms. 
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